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Abstract 



We introduce a scenario of lepton mixing in holographic composite Higgs models 
based on non-abelian discrete symmetries of the form Of = X x Zj\f, broken to 
Z2 X Z2 X Ztv in the elementary sector and to in the composite sector with 
Z^^ being the diagonal subgroup of a Z^r C X and the external Ztv- By choosing 
X = A(96) or A(384), a non-vanishing ^13 of order 0.1 is naturally obtained. We 
apply our considerations to a 5D model in warped space for the particular cases of 
X = 54,^5,A(96) and A(384) and = 3 or 5. Lepton flavour violating processes 
and electric dipole moments are well below the current bounds, with the exception 
of /X — )■ 67 that puts a very mild constraint on the parameter space of the model, for 
all presented choices of G/ . 



1 Introduction 



The data accumulated in neutrino experiments over the past years clearly show that lepton and 
quark mixing are vastly different. Several successful explanations for a lepton mixing pattern 
with two large angles and a small one in terms of a flavour symmetry can be found in the 
literature. The most prominent pattern is tri-bimaximal (TB) mixing [1] which can be elegantly 
derived with the help of the symmetries A4 [2] and 5*4 [3]. Recently, the T2K [1] and MINOS 
[3] Collaborations published indication that the lepton mixing angle ^13 is non-zero. According 
to global fits [3 El E], the best fit value of ^13 is around 0.1 0.2 and its value is different 
from zero at the 3a level [3 [7]. In the light of this, many of the models predicting TB mixing 
become disfavoured, because the deviation from ^13 = necessary to accommodate the best fit 
value of 613 0.1 0.2 is too large to be explained by sub-leading corrections. Therefore, new 
lepton mixing patterns with non-vanishing ^13 based on discrete non-abelian symmetries have 
recently been put forward The key assumptions [9j (see also [2) are that the neutrino and 
the charged lepton mass matrices are invariant under two distinct subgroups and Ge of a 
flavour group Gj, respectively, and that left-handed (LH) leptons are in an irreducible triplet 
representation of G/. Non-trivial lepton mixing is determined through the relative embedding 
of Gi, and Ge into Gf. In contrast, lepton masses remain unconstrained. 

The symmetry breaking pattern of Gf as proposed in [3l[9] is naturally realized in Composite 
Higgs Models (CHM), where the group G/ is broken to in the elementary and to Ge in the 
composite sector^ In fact, a concrete realization of this scenario was already introduced in CHM 
in [TI], for the particular case = 54 x Z3, leading to TB mixing and thus ^13 = 0. 

Aim of this paper is to generalize the scenario [TT] in order to include lepton mixing patterns 
that, like in [9], lead to non- vanishing ^13 of order 0.1 -7- 0.2. More precisely, we consider a 
set-up in which the discrete flavour group is Gj = X x Zjv, with X being a non-abelian group. 
The additional cyclic symmetry Zjv is in general needed to keep the natural explanation of the 
fermion mass hierarchy given by the Holographic CHM (HCHM) or their five-dimensional (5D) 
realizations. The pattern of flavour symmetry breaking is driven by symmetry considerations 
only, and no specific sources of flavour breaking are introduced. We focus on a scenario in which 
the Standard Model (SM) neutrinos are Majorana fermions and the type I see-saw mechanism 
explains the smallness of their masses. However, a similar analysis also applies to a scenario, 
introduced in [11], in which neutrinos are Dirac fermions. 

We choose as remnant symmetry Gi, in the elementary sector Z2 x Z2 x Zjv, while Gg of 
the composite sector is taken to be Z^^ , the diagonal subgroup of a C X and the external 
Z^r . As discussed already in [TI] , a large breaking of Gj in the composite sector is favoured for 
charged leptons, because in this way large deviations from the SM Ztt coupling are suppressed. 
At the same time, the breaking of Gj in the composite sector affecting neutrinos is required to 
be small, in order to not perturb too much the lepton mixing pattern determined by the choice 
of Gf, Ge and G^- 

^Models in warped space making use of the discrete symmetry A4, can be found in [lUj . 
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After a general presentation of our set-up in terms of 4D HCHM, we study the 5D models 
introduced in [11] in more detail for the choices (X, Z^r) = (S^jZ-^), {A5,Z^), (A(96),Z3) and 
(A(384), Z3). The full discrete symmetry of the 5D models is actually GfxY, where Y = Z'^x Z3 
is a flavour-independent factor useful to minimize the number of allowed terms. Keeping the 
prediction of the solar mixing angle 9i2 within the experimentally allowed 3a range requires that 
the flavour symmetry breaking at the IR brane for neutrinos should be smaller than 10% in all 
models, unless a Z2 exchange symmetry is imposed on the IR brane, in which case no constraint 
occurs. The corrections to the mixing angles 623 and ^13 are generically smaller. In the case 

01 X = A (96) they help to improve the accordance of the predicted and the measured value 
of 023- In the cases of X = and X = A^, in which the unperturbed value of ^13 vanishes, 
such corrections are not enough to generate a value of ^13 of order 0.1 0.2, as favoured by the 
latest experimental data and global fit analyses. Overall, the patterns derived with X = A (384) 
describe the data in the best way. Although the neutrino mass spectrum is not predicted, a 
normally ordered spectrum is preferred in the 5D models, because corrections to the solar mixing 
angle are under much better control in this case. 

For all choices of which we discuss, most of the Lepton Flavour Violating (LFV) processes 
for charged leptons are below the current experimental bounds for masses of the first Kaluza- 
Klein (KK) gauge resonances of order 3.5 TeV, roughly the lowest scale allowed by electroweak 
considerations. The main source of such processes are Boundary Kinetic Terms (BKT) for 
fermions at the UV brane. The most important bound comes from the radiative decay /i — t- 67, 
BR{ij, — )• 67) < 2.4 X 10~^^ [12], and is passed in most of the parameter space, while the expected 
future bound constrains our model. We also argue that Electric Dipole Moments (EDMs) for 
charged leptons are negligibly small. 

The structure of the paper is as follows. In section 2 we generalize the set-up of [llj to 
include generic lepton mixing patterns arising from a non-trivial breaking of a flavour symmetry 
Gf. In section 3 we apply our considerations to the 5D Majorana model of [11] and discuss the 
results for lepton mixing as well as constraints coming from charged LFV decays and lepton 
EDMs. We conclude in section 4. The relevant group theory of ^4, A^, A(96) and A(384) and 
an explicit choice of basis for their generators can be found in the appendix. 

2 General Set-up 

The set-up we consider is closely related to the one introduced in [11] for the particular choice 
of the discrete group S4 x Z3. We mainly emphasize here the key differences with respect 
to [11], referring the reader to [11] for further details. We consider in this paper only CHM 
with Majorana neutrinos, since they overall seem to perform better than the CHM with Dirac 
neutrinos, but similar considerations apply to the latter case as well. 

The Lagrangian of the CHM consists of an elementary, a composite and a mixing sector [13] : 



3 



We assume that C-m%x is invariant under the discrete flavour symmetry Gf = X x Zjy, with X 
a non-abeUan grouj which has the following features, see [S E]: i) it contains, at least, one 
(faithful) irreducible three-dimensional representation 3, as which LH leptons and right-handed 
(RH) neutrinos transform|f| and ii) it contains Z2 x Z2 and Z^r as non-commuting subgroups. The 
former requirement ensures the possibility to determine all mixing angles through the choice of 
Gf,Ge and Gi,, while the latter ensures that the resulting mixing pattern is non-trivial, because 
lepton mixing corresponds to the mismatch in the embedding of the two subgroups Z2 x Z2 and 
Zjv into X. 

The symmetry is broken in the elementary sector to G^ = Z2 x Z2 x Z^r, where Z2 x Z2 C 
X, and in the composite sector to Ge = zj^\ the diagonal subgroup of the external Zat and 
Zjy C X. Thus, all terms of C^i are invariant under G^,, while all terms of Ccomp under Ge- 

Ge is chosen as Z^^ with > 3 in order to distinguish the three generations of charged 
leptons and consequently to explain the observed hierarchy among their masses. More precisely, 
it has to be the diagonal subgroup of Z^r c X and the external Z ^ because LH leptons have to 
be assigned to a 3 of X for lepton mixing, while RH charged leptons transform trivially under 
X and carry only non-trivial charge under the external Zjy. Obviously, in order to distinguish 
among the three generations, their charges have to be different. We consider the group Z^v for 
simplicity. In a more general set-up this group can be replaced by a product of cyclic symmetries, 
such as Z2 X Z2. If the group X has three or more inequivalent one-dimensional representations, 
the additional cyclic group factor might be abandoned, because it is then possible to distinguish 
the three generations of (RH) charged leptons with the help of X alone. 

Gy consists of a Klein group Z2 x Z2 and of the external Z^v. The Klein group is the maximal 
symmetry preserved by a Majorana mass matrix in the case of three neutrinos [3], and at the 
same time it can guarantee the existence of three independent parameters corresponding to the 
neutrino masses o The external Z^v does not play any direct role for the generation of lepton 
mixing, because neither LH leptons nor RH neutrinos transform under it, but it automatically 
keeps the kinetic terms of RH charged leptons flavour diagonal in the elementary sector (see 



below) Ifl 



In the basis in which the generator Gtv of Z/^r C X is diagonal for 3 
Gn 



/cj^^ 

tj]^^ I , with ne / / , ojn = e."-"""''^ , (2.2) 



V a.]^ 



^The group X can in principle be infinite, but all examples we present in the following make use of a finite X. 

^We assume for simplicity tliat LH leptons and RH neutrinos transform in the same way under G/. 

''in the case of Dirac neutrinos Gv is not constrained to contain a Klein group, since the symmetry preserved 
by a Dirac neutrino mass matrix can be any product of cyclic symmetries which allows to distinguish among the 
three generations. The simplest case is then that Gv is a product of a cyclic symmetry Zm and the external Zjv 
with Af, iV > 3. 

^In a set-up without the external Zjv the RH charged leptons should transform as three distinct singlets under 
the subgroup Z2 x Z2 contained in Gu (and thus under X) in order to keep their kinetic terms flavour diagonal 
in the elementary sector. 
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the generators Gi and G2 of Z2 x Z2 are of the form 

Gi = ycf ''^yt^ G2 = vcf^'^v^ , (2.3) 



with 

~,diag _ n 1 n I ^diag 





^a»ag ^ 0-10, = 010, (2.4) 



and V a unitary matrix. 

The elementary sector contains three generations of SM LH and RH leptons , l'^ and three 
RH neutrinos u'^. The LH leptons and the RH neutrinos transform as (3, 1) under X x Z^r, 
while the RH charged leptons transform as {l,uj^). The elementary Lagrangian is 

Cei = limi + limR + ^Ri^^R - \{l4tM^py^^ + h.c.) , (2.5) 

where M is the most general mass matrix invariant under Z2 x Z2 x Zjv, of the form 

M = V*MdV^ , (2.6) 

with V as in (|2.3p and Md a diagonal matrix containing three independent complex parameters. 

The composite sector is an unspecified strongly coupled theory, that gives rise, among other 
states, to a composite SM Higgs field and vector-like fermion resonances ^ mixing with the SM 
fields. The fermion mixing Lagrangian Cmix is 

Cmix = -^ItK R + -r^^R^L L + -T^^R^u. L + h.C. (2.7) 

where A is the UV cut-off scale of the composite sector and and ^'"^ are fermion 

resonances transforming under G/ in the same way as and z^^, respectively. The mixing 

parameters A;^ and Xy^^ are flavour universal, because two triplets of X are coupled to each 
other, while are flavour diagonal, but non-universal, since RH charged leptons and are 
singlets under X. Integrating out the composite sector gives rise to the following charged lepton 
mass matrix (in left-right convention, ipLMiipfi): 

Mi^^p ~ ho.VHK\rM{j)''''^^''^ , (2.8) 

where vh — 250 GeV, /i is the ©(TeV) scale at which the composite theory becomes strongly 
coupled and ba are 0(1) coefficients. In this basis the charged lepton mass matrix is flavour 
diagonal and non-trivial mixing is encoded in the light neutrino mass matrix. The latter arises 
upon integrating out the RH neutrinos 
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with ba being order one coefficients. The second relation holds in the limit of universal b^. Only 
in this limit the lepton mixing matrix Upmns is given by 

UPMNS = V. (2.10) 

Obviously, a sensible choice of V implies that the resulting mixing angles are in good agreement 
with the experimental data. Deviations from the universality of 6q lead to corrections of the 
lepton mixing angles and, in order to keep the accordance with experimental data, ba generically 
have to be universal at a level of < 10% (depending slightly on the choice of V). This condition 
is equivalent to requiring that the breaking of Gf in the composite sector should be small for 
neutrinos. On the other hand, it can be of order one for charged leptons, because the values of 
ba do not have a direct impact on the lepton mixing. 

The remnant symmetry renders all couplings flavour diagonal in the composite sector 
and flavour violation is only present in (|2.5|) . All flavour changing processes are then negligibly 
small, since they are suppressed by the large Majorana mass of the RH neutrinos. The main 
source of flavour violation arises from the elementary sector, if the most general kinetic terms 
of the SM fermions compatible with the flavour symmetries are taken into account. These are 
of the form 

+ Zi)UplL + Ir{1 + Z[^)mR + ^^^(1 + Z,)ipUR (2.11) 

with Zi = VZpv\ Zy = VZj^V'^ and Zf , Zf and Z^ diagonal matrices. As explained above, 
the Ztv contained in Gu forbids flavour violating kinetic terms for IJ^. Non-trivial LFV processes 
(and further corrections to the lepton mixing (j2.10p ) are now generated and are proportional to 
the non-diagonal entries of the matrices Zi and Zy in (j2.1ip . In the limit in which the composite 
sector is G/-invariant, one can go to a basis in which the whole Lagrangian Ctot in (|2.ip is flavour 
diagonal, since in this limit all couplings and mass terms of fermion resonances \I' forming triplets 
under Gf are flavour universal. Thus, the actual amount of LFV is controlled by the size of the 
flavour violation in the elementary and the flavour non-universality in the composite sector and 
is consequently suppressed with respect to an anarchic scenario with no flavour symmetries. 

Without further constraints on the breaking of Gf in the elementary sector, the elements 
of Zf^ are expected to be uncorrelated 0(1) parameters, possibly leading to too large flavour 
violating processes and corrections to the lepton mixing. The actual effect induced by j,, 
however, depends on the degree of compositeness of and z/^ (in turn determined by the 
mass mixing terms A;^ and A,y^), since the kinetic terms of the elementary flelds always receive 
a contribution coming from the strongly coupled sector, when the fermion resonances ^ are 
integrated out (c.f. (2.7) of [IT]). Thus we have to rescale the flelds to canonically normalize 
their kinetic terms, ^ — )■ sj A^,. If the contribution from the composite sector dominates, 

S> 1, the rescaled parameters Zi^yjAi^y become suppressed. As explained in [11], Aj^^ is 
a relevant coupling and are mostly composite fields. The effect of Zy is thus negligible, 
being suppressed by the large value of A^,. The mixing term A/^, on the other hand, should be 
irrelevant, otherwise too large deviations from the SM gauge couplings of /£ would occur (see 
(|2.14p below), but it can be very close to be marginal, 7;^ ~ 0. In this case the kinetic term of 
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11 can still be dominated by the contribution coming from the composite sector due to a large 
logarithmic running from the scale A to ^, and again the rescaled parameters Zi/Ai turn out 
to be small (compare ()3.17p ). Summarizing, the flavour violating effects and the corrections to 
lepton mixing induced by Zi^i, are naturally suppressed by the dynamics in the composite sector. 

Flavour symmetries can be important also for flavour conserving observables, such as the 
EDMs. Using standard conventions, we denote the charged lepton EDMs as the coefficient of 
the dimension five operator {—i/ 2)1^(7^'^ j^l^F^i,. Being the latter a CP-odd operator, EDMs can 
be generated, if CP violation is present in the lepton sector. In absence of flavour symmetries, 
a rough one-loop estimate, assuming generic masses and complex couplings in the composite 
sector, gives 

where M^a are the charged lepton masses, Y represents a typical 0(1) Yukawa coupling involving 
the composite fermion resonances and m\j/ their typical mass. The appearance of Mi^a in (|2.12p 
is a consequence of the partial compositeness of the SM fermions which implies that the EDMs 
vanish unless both mass mixing for LH and RH charged leptons are inserted in the one-loop 
diagram. Despite the explicit dependence on the charged lepton mass, the strongest bound 
comes from the EDM of the electron, \de\ ^ 10~^^e cm [T3], giving > 101" TeV. In our 
scenario, however, the leading order result (j2.12p vanishes for canonical kinetic terms (|2.5p . 
because the relevant non-trivial phases can be removed through field redefinitionsjfl For general 
kinetic terms (|2.1ip and flavour symmetry breaking in the composite sector, the leading order 
term no longer vanishes and is estimated to be, assuming again arbitrary complex couplings in 
the composite sector, 

where 6mqi is the inter- generational mass splitting of the fermion resonances. The last two terms 
in (|2.13p are essential suppression factors that can significantly reduce the size of da- 
Let us conclude this section by showing how deviations from the SM Zlala couplings disfavour 
a small breaking of the flavour symmetry in the composite sector for charged leptons, 6q <C 1 
in (|2.8p . For simplicity we set Zi = Zf = Z^ = 0, their effect being sub-leading. The coupling 
deviations arise from mixing of the SM leptons with fermion resonances, induced by the mixing 
terms in (|2.7p . One schematically has 



9lc^ ^ rnl'"'^' \AJ 
Using the charged lepton mass formula ()2.8p . we can write 

^9l^L ^9l^R ^la 1 



(KrU) . ^ = L,R. (2.14) 



9i^L 9i^R ml ml hi ' 



(2.15) 



^This is not true in general, if one considers non-minimal scenarios, in which for example SM fermions mix 
with more than one state of the composite sector. 
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The tension between having parametricahy smah ba and sufficiently small gauge coupling devi- 
ations is obvious from ()2.15p . This is particularly important for the tau lepton due to its larger 
mass. Deviations from the SM gauge couplings for (charged and neutral) leptons have been 
constrained by LEP at the per mille level [15] • For fermion resonances at the TeV scale, we see 
that the right-hand side of ()2.15p is below 10~^ for b^- slightly below one. Notice that one might 
actually cancel the leading term ()2.14p in 6gi^j^ for either charged leptons or neutrinos using 
appropriate symmetries [TB], but not both at the same time. The relation (|2.15p is then always 
valid for at least one of the two components of the LH doublet 

3 5D Realizations 

Models with gauge-Higgs unification in warped space based on an SO{5) x U{l)x gauge symme- 
try [13l [T7] are an explicit and particularly interesting weakly coupled description of the scenario 
outlined in section f2V\ We consider fiavour groups of the form: 

GfXY. (3.1) 

The additional discrete symmetry Y, 

Y = Z',xZ'i, (3.2) 

is introduced in order to minimize the number of couplings in the bulk and at the branes. We 
present four examples in the following leading to different results for the lepton mixing angles. 
Our first choice of G/ is 

G/ = 54 X Z3 (3.3) 

which is discussed in [11] and leads to TB mixing, see (j3.10p . We repeat its analysis here for 
completeness. The choice 

Gf = ^5 X Z5 (3.4) 

gives rise to the so-called Golden Ratio (GR) mixing in which the solar mixing angle 612 is 
determined in terms of = ^(1 + -v/5) |19] . while vanishing 613 and maximal 623 are predicted, 
see (|3.1ip . We include its discussion, albeit it might seem to be disfavoured due to the prediction 
613 = 0, in order to show an example in which the external Z^r factor is different from Z3. As 
third and forth choice, we discuss the cases 

Of = A(96) X Z3 , Gf = A(384) x Z3 , (3.5) 

since it has recently been shown that these groups can naturally lead to ^13 ~ 0.1 0.2 [9]. 
They give rise to two inequivalent mixing patterns each (differing in the value of the angle ^23); 
and thus we get in total four different possibilities. Following the notation of [9], we call them 
Ml and M2 for A(96), see (fXT2|) and ([XT3]1 . and M3 and M4 for A(384), see (l3Tl]l and (i3T5]l . 
respectively. 

^Notice that warping is not a necessary ingredient. For viable models in flat space see [18) . 



8 





Bulk 


UV 


IR 


GfxY 


Gf,\jv X Yuv 






(3, l,UJs,^3) 


(1,-1,1,^^3) 
(-1,1,1,^3) 
(-l,-l,l,a;3) 










«",'^3) 




(3,1,W3,1) 


(1,-1,1,1) 
(-1,1,1,1) 
(-1,-1,1,1) 


(^]^",W3) 



Table 1: Transformation properties of the 5D multiplets (,e,a and ^^^a under Gj x Y and 
the subgroups Gj^uv x Yuy and Gjjr x Iir. The values of and for each group can be 
found in table [2j 

The flavour symmetry is broken at the UV and IR branes to 

G/,uv X Yvv = Z2 X Z2 X Zjv X Z'3' , G/,ir x Y^ = Z^^^ x Z'3 , (3.6) 

with iV = 3 in the case of X = S'4, A (96) and A (384), and iV = 5 for A5. 

The lepton particle content of the model is identical to the one of [11] with respect to the 
gauge group: three 5D bulk fermions ^/^a, ^e,a and £,u,a, in the fundamental, adjoint and singlet 
representations of 50(5) are introduced (for details of the notation see [H]). All of them have 
vanishing U{l)x charge. Their flavour properties are reported in tables [T] and [21 

The most general Gjjr x Yi^ invariant mass terms at the IR brane are 

— j ^ (jnm,a {Ll,aLL2,aR + LahLaRj + "1^ „ I'aL'^aR + h.C^ . (3.7) 

^ a=e,fj,,T 

In the particular bases chosen for the different groups G / , see section [2] and the appendix, these 
terms are flavour diagonal. The fields La and Oa are components of the 5D multiplet S^i^a, 
L2,a and La are contained in ^e,a and Va = ^u,a- The only G/^uv x Yuv invariant mass terms 
at the UV brane are Majorana mass terms for RH neutrinos: 

- >CuV = ^ T^aR-^VV,apVpR + h.C. (3.8) 

with 

Mvw = V^muYV^ , (3.9) 
muv = diag (muv,e , ^^\]V,^l■, "^uv.t) and V as in ()2.3p . Up to (removable) phases and signs, 
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X 


Sa 


^5 


A(96), Ml 


A(96), M2 


A(384), M3 


A (384), M4 


N 


3 


5 


3 


3 


3 


3 


ria 


(0,2,1) 


(0,1,4) 


(2,1,0) 


(2,0,1) 


(1,2,0) 


(1,0,2) 



Table 2: Values of N and for the different choices of non-abelian discrete groups X and 
mixing patterns. 



the explicit form of V is as follows: 

/ 

X = S4 '. V = Utb — 



COS 
sin e^2^ 

sin 6*^2^ 
\ V2 



- sin 6f^ 

cos 6*^2^ 



V2 



COS 6*^2^ 



V2 



^^ 

1 



with t, 



,an 



(3.10) 



(3.11) 



X = A(96), Ml : V 



1 



/_i(^+l)l i(V3-l) 
i(^/3-l) 1 -i(^/3 + l) 



1 1 1 



(3.12) 



X = A (96), M2 : V equal to the one of Ml with 2^"^ and 3"'^ rows exchanged , (3.13) 



X = A(384), M3 : y = ^ 

v3 



/-iV4 + V2 + V6 1 -iV4- \/2- V6\ 
i\/4T72^V6 1 -i\/4- ^/2 + 



V 



1 



(3.14) 



X = A(384), M4 : V equal to the one of M3 with 2"^ and 3'''^ rows exchanged . (3.15) 

In order to discuss the result for lepton mixing analytically, we consider charged lepton and 
neutrino mass matrices in the Zero Mode Approximation (ZMA), including the effect of the 
dominant flavour violating BKT 

JObkt = Ll{x, R){RZi)^Ll{x, R) , (3.16) 

with Zi being constrained by Z2 x Z2 x Z^v to be of the form Zi = V(iia.g [z^u z^,hZrl)V^% The 
effective BKT relevant for the single KK modes are obtained by multiplying Zi with the square 



*The coefficients Zei, z^i and Zt-i were denoted by Zei, z^i and Zri in [TT]. We use the latter notation for the 
entries of Zi defined in p.l7p . 
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of their wave function profile evaluated at the UV brane. For zero modes we get, if the bulk 
mass parameter q of the fermions ^^^q fulfills q = 1/2 + (5c, at linear order in 5c, 

Zi ^ (log-i ^ + Sc)Zi ^ + Sc)Zi , (3.17) 

taking R of the order of the inverse of the reduced Planck mass and R' of the order of the inverse 
of the TeV scale. The matrix Zi in ()3.17p should roughly be identified with Zi introduced in 
()2.1ip . and 1/35 is the 5D counterpart of the suppression factor coming from the composite 
sector, discussed in section 2. The latter plays a crucial role in naturally suppressing most of 
the LFV processes well below their current experimental bounds and in keeping the corrections 
to the lepton mixing small. The charged lepton mass matrix in the ZMA reads, after canonical 
normalization of the kinetic terms (j3.16p and additionally rotating LH charged leptons with 
(again all relevant notation can be found in [llj ) 

Ml = fcj^c, I V , ^ V^"^ I . (3.18) 



a/3 

The light neutrino mass matrix, after integrating out the heavy RH neutrinos and canonically 
normalizing the kinetic terms (this time without rotating LH neutrinos with V"^), IS 

2i?'2 i^Rj \^1 + ZP ^/P "^UV VP ^l + Z^I ^ 
In this basis the charged current is of the form 

IlW-Vvl (3.20) 

which coincides with the result given in ()2.10p in the limit in which the BKT are set to zero 
and the mass parameters rn^^ ^ as well as the factors pa are taken to be universal. In the 
phenomenologically interesting region of the parameter space in which the bulk mass parameter 
Q of ^i^a is close to 1/2, the mass parameters mj^a '^^'^ order one without affecting 

considerably the universality of the parameters pa and thus the results for lepton mixing. On 
the other hand, we still need to assume a small breaking of G/ at the IR brane in the neutrino 
sector, i.e. 

KR,a = KR,o(l+'5a)> (3.21) 

with \6a\ <^ 1, in order to keep their impact on the mixing angles under control. As explained 
in [11] , the parameters in ()2.9p should be identified with the mass parameters mj^p^ ^ and thus 
6a measure the non-universality of 6q,. As we see below, \5a\ ^ 0.1 are required, reflecting that 
the parameters 5q, have to be nearly universal. We analyze this issue in more detail in the next 
subsection. Alternatively, we can require the invariance of the IR localized Lagrangian under a 
Z2 exchange symmetry, under which 

l'a{x,R') ^ Ua{x,R') , Li^a{x,R') ^ L2,a{x,R') , La{x , R') ^ La{x , R') , (3.22) 
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so that we can take w-n^^ = "ijj^^ = 1 (up to an irrelevant sign per generation). We denote 
this constrained model as Z2-invariant model. 

3.1 Lepton Mixing 

We first discuss the phenomenological constraints on the size of (5q,, as defined in ()3.2ip . in the 
ZMA at linear order in the perturbation Sa- In doing so we set Zf = and neglect the non- 
universality of the parameters pa, which is small for q ~ 1/2. Then we analyze these constraints 
numerically by taking into account the first KK level, still without considering the effect of BKT. 

We find the following analytical results, for normally ordered neutrinos with a lightest neu- 
trino mass rriQ = 0.01 eV and solar and atmospheric mass square differences Ajn^^j = 7.59 x 10~^ 
eV^ and Aml^^ = 2 AO x 10~^ eV^ [20]: in the case of ^4 and A^, deviations from ^13 = and 
maximal 823 are proportional to the breaking of /.t — r symmetry (5^ — 6r) and are nearly the 
same (the values in square brackets, if given, refer to A5): 



sin 013 0.05 — 6r \ , 
1 
2 



sin^ 023 -1 + 0-82 [0.83] {6^ - 6r) , (3.23) 



sin^ 6*12 w ^ [0.28] + 1.58 [1.43] (2^ - 
For the mixing pattern Ml arising from A (96) we find 

sin^ « 0.04 + 0.135e - 0.115^ - 0.035^ , 

sin^ 6*23 ^ 0.65 - 0.025e + 0.765/, - 0.74 6r , (3.24) 
sin^ 6^12 w 0.35 + 3m5e - 0.675^ - 2.425^ . 

The results for M2 are related to these by exchanging 6^ and 5r and by replacing the unperturbed 
value of sin^ 023 by 0.35 together with a sign change in its corrections (remember that Ml and 
M2 are related by the exchange of the second and third rows of V). As can be seen, for 5^ ~ 6r 
corrections to the atmospheric mixing angle become suppressed. For the mixing pattern M3 
coming from A (384) we get 

sin^ 6*13 « 0.01 + 0M6e - 0.0045^ - 0.035^ , 

sin^ 6*23 ~ 0.42 + O.OUe + 0.805^ - 0.81 6r , (3.25) 
sin^ 6*12 ^ 0.34 + 3.155e - 2.025^ - 1.135^ . 

The mixing angles and their corrections in (5q- associated with the patterns M3 and M4 are related 
in the same way as those of Ml and M2 (the unperturbed value of sin^ 023 is 0.58 in the case of M4). 
Again, 023 only receives small corrections for 6^ ^ 6r- As ()3.23p - ()3.25p show, the corrections 
to the solar mixing angle, the one which is experimentally determined with best precision, are 
generally the largest with the coefficients of 5a being larger than one. This also implies that the 
above perturbative expansion makes sense only for \6a\ ^0.1. As mentioned in |11] . the validity 
of the expansion in 5a strongly depends on niQ and gets worse for increasing uiq. In the case 
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-0.10 -0.05 0.00 0.05 0.10 

Figure 1: The solar mixing angle 6*12 as function of the deviation S, parametrizing the non- universality of 
the masses mj'j^ ^ , see p.26p . The horizontal gray band corresponds to the 3cr range as given in [S] using 
new reactor fluxes. The different symbols and colors distinguish the various models: 6*4 (red points), 
(blue crosses), A(96) and mixing pattern Ml (light green full triangles), A(96) and mixing pattern M2 
(dark green open triangles) , A(384) and mixing pattern M3 (black full diamonds) and A(384) and mixing 
pattern M4 (violet open diamonds). The curves for Ml and M2 and M3 and M4, respectively, lie on top of 
each other due to the parametrization chosen in p.26p . The plots refer to q = 0.52, = —0.365, h = 1/3 
and normal neutrino mass hierarchy. The masses muv.a a-re chosen such that the lightest neutrino mass 
is Too = 0.01 eV and the values of the solar and atmospheric mass square differences Ato^^j = 7.59 x 10^^ 
eV^ and Ato^^^ = 2.40 x 10"'^ eV^ [20] are reproduced using the ZMA. For simplicity, we take the IR 
mass terms toJjj „ = 1 and set all BKT to zero. We vary m'^^ q, see p.2ip . between 0.3 and 1. 



of inverted mass hierarchy the coefficients multiplying the linear perturbations in 612 are more 
than one order of magnitude bigger than the ones in (j3.23p - p.25p . implying that a perturbative 
expansion in 5a is not valid for any value of tuq. This behaviour is in general expected due to the 
near degeneracy of the two heavier neutrinos in the case of an inversely ordered mass spectrum. 
Note that in the limit of universal 5a all corrections to the mixing angles vanish. 

In our numerical analysis we discuss the maximal allowed size of the corrections 5a in order 
to keep accordance with experimental data. We choose a particular parametrization for the 
deviations of mj'j^ ^ from universality in terms of only one parameter 5: 





54, ^5 


6e 


= 0,<5^ 


= S,5r 


= 


A(96),M1 


A(384),M3 


Se 


= S,5^ 


= (5, 5r 


= 


A(96),M2 


A(384),M4 


Se 


= S,5^ 


= 0,5r 


= 5 



This particular parametrization leads to the same results for the mixing angles ^13 and 612 in 
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case of the patterns Ml and M2 as well as for M3 and M4. The atmospheric mixing angle acquires at 
the same time a correction which is the same in size, but opposite in sign for Ml (M3) and M2 (M4). 
As already obvious from the analytical results, the corrections are the largest for the angle 9i2 
whose dependence on 6 we report in figure [U together with its experimentally allowed 3a range 
using the new estimate of reactor anti- neutrino fluxes [6]. Clearly, \5\ has to be smaller than 
0.07 in all cases. Depending on the unperturbed value of sin^ 9i2, either positive or negative 5 is 
better compatible with the data for increasing \5\. All cases apart from X = 84^ prefer 6 <0 for 
the particular choice (j3.'26p . We do not show our results for the other two mixing angles, since in 
the case of ^13 all corrections for \6\ ^0.1 are small: in the case of 5*4 and sin^ ^13 ^ 4 x 10~^ 
holds, compare (|3.23p . while for patterns Ml and M2 we find 0.041 < sin^ 6*13 < 0.047 and for M3 
and M4 we get 0.0084 < sin^ ^13 < 0.015. We clearly see that corrections associated with the non- 
universality of the neutrino Dirac mass terms are not sufficient in the cases X = S4 and X = 
to explain ^13 ~ 0.1-^0.2. For \5\ < 0.1, the corrections to the atmospheric mixing angle always 
keep sin^ 6*23 in its 3a range [6] for S4, A^ and A(384). Since the unperturbed value of sin^ ^23 
is at the edge of the 3a range for the patterns Ml, sin^ ^23 ~ 0.65, and M2, sin^ ^23 ~ 0.35, 
corrections with 6 < for the parametrization ()3.26p are welcome because they improve the 
agreement with the results from global fits. At the same time 612, see figured! remains in its 
experimentally allowed 3a range for negative 5 with \6\ < 0.04. Generally speaking, the patterns 
MS and M4 are the most promising ones even taking into account corrections coming from the 
non-universality of the neutrino Dirac mass termso In the case of the patterns Ml and M2 these 
corrections help to improve the accordance with the experimental data; however, these patterns 
are not favoured by the latter. 5*4 and A^ mainly fail to give a good fit to the data because of 
the too small value of ^13. 

The effects of BKT and of non-universal p on the lepton mixing, neglected in the above study 
being sub-leading, become relevant in the Z2-invariant model, where 5a vanish. In the latter 
model we find numerically for S4 and for A^ (see caption of figure [2] for details on the chosen 
parameters) 

0.32 < sin^ 6^12 < 0.35 and 0.27 < sin^ 6^12 < 0.29 , (3.27) 

respectively, and 

0.48 < sin2 023 < 0.52 , (3.28) 

showing that sin^ ^12,23 only get corrected by less than 0.03. Corrections to ^13 are negligible 
and thus 6*13 ~ 0.1 0.2 [3 El [H] cannot be achieved in models with S4 and A^, by taking into 
account the BKT or the deviation of the parameters p from universality. For A(96), Ml [M2] we 
get 

0.63(0.32] < sin^ 6^23 < 0.69(0.36] , (3.29) 

^There is a slight dependence of which of the two patterns performs best on the used global fit analysis: using 
[6l[8] we find pattern M3 to be the best one, while [7] prefers M4 over M3. This difference originates from the fact 
that in [3 the best fit value of sin'^ 623 is larger than 0.5, whereas it is smaller than 0.5 in 
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Figure 2: Branching ratio of yu — > 67 as a function of the UV BKT (5z = 3{Zi)efi in the Z2-invariant 
model. The continuous and dashed lines are the current {BR{ii — > 67) < 2.4 x 10^^^) and the expected 
future bound {BR{fj, ej) < 10^^'^) given by the MEG Collaboration [12], respectively. The parameters 
ci, Cu, h, niQ and the neutrino mass hierarchy as well as the color coding are chosen as in figure [TJ Note 
that no constraints coming from lepton mixing are taken into account in this data set; however, the 
constraints coming from gauge coupling deviations are satisfied by each point plotted. 

keeping the atmospheric mixing angle compatible with data only at the 3a level [6], as well as 
0.34 < sin^ 6*12 < 0.37 and 0.04 < sin^ 013 < 0.053[0.047] . (3.30) 
For A (384), M3 [M4] we find analogously 

0.40(0.56] < sin^ ^23 < 0.44[0.59] , (3.31) 

respectively, and 



0.32 < sin^ 012 < 0.35 and 0.01 < sin^ 6*13 < 0.012 . 



(3.32) 



3.2 LEV Processes 



In this subsection we examine the bounds coming from LFV processes in all six models proposed 
(see ^T] for early analysis of LFV bounds in warped models). As discussed in detail in |llj and 
repeated in section 2, these processes are non-negligible only when the leading flavour violating 
BKT ()3.16p at the UV brane are taken into account. The tree-level decay fi — t- 3e, fi—e conversion 
in nuclei and — )• 67 depend quadratically on the off-diagonal entry (efi) of the matrix (j3.17p . 
The relevant combination of Zai varies from case to case because Zi depend on the mixing matrix 
V. We plot in figure [2] the bounds arising from ^ — )• 67 as a function of 5z = 3{Zi)e^ for the 
Z2-invariant model0 obtained by a numerical computation in which the first KK mode of each 

^'^The factor 3 in the definition naturally arises in the ^4 case, and is left to match the convention used in [11] , 
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tower of states has been kept. BR{^ — )• 67) is essentially the same in all models and shows 
the expected quadratic dependence on 5z. The branching ratio is always below the current 
experimental limit for < 1 and the expected future MEG bound BR{^ — )• 67) < 10~^^ will 
require \5z\ < 0.2. Similar results hold for the model with no Z2 exchange symmetry since the 
IR masses ^ do not play an important role. Finally we note that in all models the branching 
ratios of /x — 3e, /x — e conversion in Ti as well as of the radiative r lepton decays are below 
their experimental bounds for \5z\ < 1. 

3.3 EDMs 

Lepton EDMs in the 5D model are completely calculable, because gauge invariance forbids the 
appearance of uncalculable 5D bulk or boundary operators that reduce, upon KK reduction 
in 4D, to the EDM operators [—i /'i)la(y^^ ^^laFf^u ■ The size of da is negligible, due to the 
relatively few sources of complex parameters. The bulk parameters q, Cq and c^, are real, while 
the IR and UV mass terms contain in total nine complex phases. Those of the IR mass terms 
'^iRo ~ l"^mal^*^" be removed by a simple field redefinition: 



Note that the terms in the bulk Lagrangian do not change if the transformations (|3.33p are 
applied, because the terms involving the 5D fields ^i^a and ^^^a are invariant under the flavour 
symmetry t/(3) and those containing £^(,,a under ^7(1)"^. In the field basis of (j3.33p . all phases are 
encoded in the UV Majorana mass terms: the phases associated with the UV mass terms muv,a 
as well as the phases 9^. As has been shown in [TT] (cf. (4.34)) by considering the KK expansion 
of the RH neutrinos VaR, the UV localized Majorana mass term gives effectively rise to only 
one heavy RH neutrino per generation, with mass of order 10^^~^^ GeV, while the remaining 
orthonormal combinations of KK states are not sensitive to the UV Majorana mass term and 
have masses setting in at a few TeV. As consequence, the EDMs induced by the UV Majorana 
mass terms are mediated by these three heavy states and are completely negligible. Notice 
that the transformations ()3.33p do not involve the fields ^i^a and hence no phases appear in the 
BKT ([3T6]) . The next-to-leading flavour violating BKT at the UV brane is Zy in ([2TT]) . which 
acquire non-trivial phases after the transformations ()3.33p . Due to the field localization of 
in the extra dimension, however, the effective BKT are strongly suppressed, Zy < 10~^, if values 
of Cy are used which are suitable for reproducing correctly the scale of light neutrino masses, 
e.g. Cy = —0.365 like in figures [T] and [2] (this is equivalent to the suppression mechanism of Zy 
explained below (|2.1ip for a composite z^|j). Using the estimate (|2.13p with Zi being replaced 
by Zy, we find that lepton EDMs are well below the current experimental bounds. 

4 Conclusions 

We have extended the class of 4D HCHM based on the non-abelian flavour group 5*4 X Z3 
introduced in [11] to general non-abelian discrete groups of the form X x Z^r- This allows to 





(3.33) 
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consider HCHM which predict a promising lepton mixing pattern with non-vanishing ^13 ~ 
0.1 0.2 [9] as favoured by current data. In a 5D scenario with Majorana neutrinos, we have 
computed in detail the lepton mixing for four particular choices of X x Z^y. We have shown that 
flavour symmetry breaking effects at the IR brane affecting neutrinos can be at most 10% in 
order to not perturb too much the original predictions for the mixing angles. We have argued, 
like in |llj . that a Z2 exchange symmetry can be imposed on the IR brane avoiding these 
restrictions. All LFV processes and EDMs are below the experimental bounds. 

As discussed in [llj . the typical mass scale of the vector- like leptonic fermion resonances in 
the 5D models is around 2 TeV. The prospects to produce and observe these resonances at the 
LHC are unfortunately quite limited. A possible signature might be the observation of the decay 
^ — 7- 67, considering that its typical branching ratio is within the range of the expected future 
bound of the MEG experiment. 

It would be interesting to extend our considerations to the quark sector and discuss quark 
mixing in a similar manner. 

Acknowledgment s 

M.S. would like to thank Andrea Romanino for useful discussions and Kaustubh Agashe for a 
useful e-mail correspondence. C.H. would like to thank Ferruccio Feruglio for useful discussions 
and the Aspen Center for Physics for kind hospitality during the preparation of this work. 

A Group Theory of ^4, A5, A(96) and A(384) 

In this appendix we report some details of the group theory of S'4, A^, A (96) and A (384). All 
symmetries can be defined in terms of two generators S and T. We show the relations the latter 
have to fulfill in order to generate one of the groups, S'4, A^, A (96) and A (384), and give their 
explicit realization for an irreducible triplet. The basis chosen for S and T is such that the 
Lagrangians introduced in (j2.5p . (j3.7p and (|3.8p are reproduced. In the following we report the 
explicit form of S. The generator T can be easily computed using ()2.2p and table [21 once the 
product of S and T giving Gn is known. 

In the case of S'4, S and T fulfill the relations |22] 



52 = 1 , T^ = l , {STf = 1 



(A.l) 



The explicit form of S can be chosen as 



1 



/ 1 -2 -2\ 



S 



3 



-2 -2 1 



(A.2) 



V-2 1 -2/ 



The generators of Ge and Gi, are the following 



GN = iSTf, Gi=S, G2 = {ST^f. 



(A.3) 
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The explicit form of T is then 



T = SG^ 



N 



1 

3 



/ 1 -2oj-i -2ojI 
-2 -2uJs 



(A.4) 



\-2 6^3 



-2ioi 



We also briefly comment on the relation between the generators S and T used here and the set 
of generators S, T and U used in [11] in order to describe the group S^. One can check that T 
is similar to Gn = (STf, 5 to G2 = {ST'^f and U to G1G2 = T'^ST^. 

The group is generated through S and T being subject to the conditions [TS] 



We choose the realization of S to be 



1 



{STf 



1 . 



(A.5) 




(A.6) 



with = (1 + \/5)/2 for one of the irreducible triplets. In this case the remnant subgroups Ge 
and Gi, are generated through 



Gn = T, Gi = S, G2 = T^ST^ST^ . 
The group A (96) is generated through 5 and T fulfilling the relations [9] 
52 = 1 , T^ = l , {STf = 1 , {ST^^STf = 1 . 



(A.7) 



(Ai 



In order to realize the mixing pattern called Ml, the most convenient choice of basis for S for a 
faithful irreducible triplet is 



/-I + VS -1 -1-^3^ 

-1 -1- V3 + 
(^-1-^/3-1 + ^/3 -1 



(A.9) 



For the other mixing pattern M2 it is convenient to choose S like in ()A.9p . however with second 
and third rows and columns exchanged, respectively. This is clear because the mixing matrices 
V associated with Ml and M2 are related by the exchange of the second and third rows. In- 
dependently of the mixing pattern, the remnant subgroups Ge and G^ are generated through 

m 

Gn = ST, Gi = S, G2 = {ST^f. (A.IO) 
Similarly, the group A (384) is generated with S and T fulfilling [9] 



1 



il6 



1 , {STf = 1 , {ST-^STf = 1 . 



(A.U) 
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Again, it is convenient to choose two different bases for a faithful irreducible triplet representation 
in order to generate the mixing patterns M3 and M4, respectively. For a model incorporating the 
pattern M3 we choose 



6 



/-2 + \/2 + ^/6 -2(1 + ^/2) -2 + ^/2-^/6 

-2(l + \/2) -2 + \/2-\/6 2(-l + \/2 + ^/3) | . (A.12) 
V -2 + ^2 - \/6 2(-l + v/2 + V3) -2(1 + ^2) 



For an explicit model leading to the pattern M4, the basis in which S is like in ()A.12p with second 
and third rows and columns exchanged, respectively, is the most appropriate one. Again, this is 
obvious considering the relation of the mixing patterns M3 and M4. Independently of the mixing 
pattern, the remnant subgroups Ge and are generated through [U] 

Gn = ST, Gi=S, G2 = {ST^f. (A.13) 
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